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Abstract 

We analyze the diffraction regime of the Kapitza-Dirac effect for par¬ 
ticles entangled in momentum. The detection patterns show two-particle 
interferences. In the single-mode case we identify a discontinuity in the set 
of joint detection probabilities, associated with the disconnected character 
of the space of non-separable states. For Gaussian multi-mode states we 
derive the diffraction patterns, providing an example of the dependence 
of the light-matter interaction on entanglement. When the particles are 
identical, we can explore the relation between exchange and entangle¬ 
ment effects. We find a complementary behavior between overlapping 
and Schmidt’s number. In particular, symmetric entanglement can cancel 
the exchange effects. 

Keywords: Diffraction by light gratings; Two-particle interference; Entan¬ 
gled states; Exchange effects; Discontinuous processes 

1 Introduction 

Entanglement leads to physical effects unattainable for product states. One of 
the best known examples of such effects is multi-particle interference, which has 
been extensively studied in both the theoretical and experimental realms [1, 2]. 
In general these studies are within the framework of quantum optics. It would 
be interesting to consider similar schemes for massive particles. One natural 
candidate is the interaction of massive particles with light gratings. In particu¬ 
lar, the Kapitza-Dirac effect, proposed long ago [3, 4], has been experimentally 
verified for atoms [5] and electrons [6]. 

We shall consider in this paper two-particle diffraction by light gratings 
with the particles entangled in momentum, showing the presence of two-particle 
interference effects absent for factorizable states. The arrangement reflects how 
the light-matter interaction is modified by entanglement. Other examples of the 
dependence of the interaction on entanglement can be found in the literature 
[7, 8], 
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Similar interference effects have been already discussed for the other type 
of non-product states, the (anti)synnnetrized states of identical particles [9]. In 
the case of single-mode states our work bears many mathematical resemblances 
to that paper, but here we shall consider some aspects of the problem not 
discussed there. We shall show the existence of a discontinuity in the set of 
joint detection distributions when the two momenta are equal. This property 
reflects the discontinuous nature of the space of entangled states. 

For multi-mode states our approach departs from the elementary qualitative 
treatment given in [9]. Here, we consider entangled Gaussian states and we 
can derive exact analytical expressions for the joint detection patterns. These 
patterns, corresponding to the sum of a large number of Gaussian terms, can 
accurately be described by a single Gaussian. 

The main physical novelty in the multi-mode case occurs for identical par¬ 
ticles. Then we can observe at work three quantum processes (diffraction and 
exchange and entanglement effects) at once. Our first finding in this scenario 
is the existence of a type of complementarity between overlapping (the variable 
ruling the intensity of the exchange interaction) and entanglement degree. When 
one of them increases the other must decrease. To evaluate the entanglement 
degree we use the Schmidt number [10, 11], which is a measure well-suited for 
continuous variable problems. Another consequence of the physical connection 
between particle identity and entanglement is the cancellation of the exchange 
effects for symmetric entangled states. 


2 Single-mode states 

We consider in this section single-mode states. First we describe the arrange¬ 
ment and the equations describing the light-matter interaction. In the second 
subsection we derive the interference patterns in the position representation. Fi¬ 
nally, we discuss the singular behavior of the system when the initial momenta 
are equal. 


2.1 The arrangement and fundamental equations 

A sketch of the arrangement can be seen in Fig. (1). A source generates pairs 
of particles in entangled states. Each particle interacts with a light grating. 
The gratings, two standing light waves formed by counter-propagating lasers 
with different wavelengths, are denoted as L and R. After the gratings we place 
detectors working in coincidence to determine the joint probabilities. 

The only relevant variables in the problem are those parallel to the grating 
[4], reducing the description of the system to a two-dimensional one (one variable 
for each particle). The initial entangled state in momentum is 


\<Po >= 


V2 


(b >l \q >r +\q >l \p >r) 


(i) 
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where p and q denote the longitudinal momenta of the particles and the sub¬ 
scripts L and R refer to the grating with which they interact. Next, we must 
describe the interaction of the particles with the gratings. The process is ruled 
by the liglitshift potential V = Vq cos 2 Kx, with K the wavenumber of the laser 
light and x the spatial coordinate in the direction parallel to the grating [4]. 
We assume that the potential intensity is the same in both sides but the laser 
wavenumbers, K R and K R , differ. In the diffraction regime of the Kapitza-Dirac 
effect the Raman-Nath approximation, neglecting the kinetic part of the mo¬ 
tion, holds [4]. Then, as the initial wave function for the mode p in the position 
representation is i/jq = e tpx ^ h , the wave function after the interaction will be 
tp T = e~ lVT / h ipo, with r the interaction time. The exponential can be rewritten 
using the expression e zzcosv = J n (z)e lnv> , with J n the Bessel function of or¬ 
der n. Using simple trigonometric relations we obtain ip T = b n e‘ ( ' 2nhK+p ^ x ^ h 
with b n = i n e~ lw J n (—w) and w = Vqt/2Tl. 



Figure 1: Arrangement for two-particle Kapitza-Dirac diffraction. The source 
S produces pairs of entangled particles that interact with the light gratings L 
and R. The continuous, dashed and dotted lines after the gratings correspond 
respectively to particles with n = 0, ±1 and ±2. 


Returning to the momentum picture this evolution can be expressed as 
|p >l~> J2n bn\p + 2nfi.A'i >l- The final state of the complete system is 


\4> >= — 7 = ^2 b n b m (\p + 2nhK L > L \q + 2mTiK R > R 

' nrt rr» — — rv"s 

+ |g + 2 nhK L > L | p + 2mhK R > R ) 


n,m =—oo 


( 2 ) 


From this expression we can derive the joint detection probabilities. The prob- 
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ability of measuring the system in the state |p + 2 tiUKl >l \q + 2mhKR >r 
is 

P(rip,mq) = | L <p + 2nhK L \ R <q + 2m‘hK R \<j) > | 2 = i|6„6 m (5(0)5(0) + 

oo 

b n ib m i5(p + 2nhK R — q — 2n'hKL)S(q + 2mUK R — p — 2m'hK R )\ 2 (3) 

n / ,m / =—oo 

where we have used the orthogonality of the momentum states < p|p* >= 
d(p — p*)- The above expression must be interpreted in terms of distributions 
where, for instance, b n S(0) = b n < p + 2 uTiKl\p + 2nTiK R >. When the two 
conditions 

2 (n — n')hK R = q — p (4) 

and 

2 (m — m')hK R = p — q (5) 

simultaneously hold for some values of n, m, n! and m! we have interference 
effects. In effect, when these conditions are not fulfilled, (n — v!)K R ^ ( m’ — 
m)K Rl the detection probability is P(?Zp,m^) = \b n b m \ 2 /2 with no interference 
effects. In contrast, when the conditions hold the probability becomes 

P(rip,mq) = ^\b n b m \ 2 + ^\b n >b m >\ 2 + Re{b* n b* m b n >b m ') ( 6 ) 

The third term on the r. h. s. corresponds to the interference term between the 
alternatives (i) n photon absorptions in L and m in R, and (ii) n' absorptions in 
L and m' in R. The indistinguishability of both alternatives implies that their 
probability amplitudes must add, giving rise to interference effects. 

We note that, for given K R and I\ R , only in some cases we have interference 
effects. The two above conditions lead to the relation (n — n')K L = (m r ~m)K R 
that can be only fulfilled when K R /K R is a rational number as, for instance, in 
the case K R = K R . On the other hand, for each pair of values p and q (with 
the exception of, see later, p = q) there are always pairs of laser wavevectors 
showing interference. For instance, taking n = m = 0, n' = — 1 and m! = 2 
there are interference effects for K R = (q — p)/2Ti and K R = 2K R . 

2.2 Interference in the position representation 

In the position representation the initial entangled state is d'o = + 

' i /’i(2/)V’fl( a: ))/V / 2, with x and y the longitudinal spatial variables of the two 
particles. After the interaction this expression becomes 

*r(®, y) = -^(Mxy px,h v R {yV qv/h + MyV qy/n Mxy px/h ) (7) 
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with = ^2 n b n exp(i2nKLx) and a similar expression for tpn. The final 

two-particle detection probability, |4 r r (a;, y)\ 2 , is 

P(x, y) = P pro (x, y) + Re(ipl(x)^ R {x)^(y)ip L (y)e l{q - p){x - y)/n ) (8) 

with 2 P pro = \pL(x)\ 2 \tp R {]j)\ 2 + \ip L (y)\ 2 \<p R (x)\ 2 the probability of the product 
state corresponding to a mixture of the initial product states ipL{x)i/jR(y) and 
^Liy^Rix) with equal weights 1/2. 

The second term on the r. h. s. of this equation represents the two-particle 
interference effects, with the typical dependence on trigonometric functions of 
x — y. At variance with momentum interference there are no constraints be¬ 
tween the values of p,q,K R ,KR and the n's. As a matter of fact, the spatial 
interference contains all the possible numbers of photon interchanges (all the 
n's) through Lp. Moreover, there is spatial interference for any value of the ra¬ 
tio KijK r. The above properties show the different qualitative behavior of 
interference in the position and momentum representations for the same initial 
entangled states. 

The interference effects discussed in this subsection and the previous one are 
mathematically similar to those in [9] for identical particles. In that paper they 
were considered in detail and we shall not continue here that line of argumen¬ 
tation, which can be easily translated to our problem. Instead, we shall analyze 
another physical aspect of the problem that was not treated there, the existence 
of a discontinuity. 

2.3 Discontinuity 

Let us consider the particular case p = q in Eqs. (4) and (5), the conditions for 
the existence of interference effects. In this case the effects only exist for n = n' 
and m = ml. However, the resultant patterns cannot be considered as genuine 
two-particle interference as the necessary condition to have interference is the 
existence of indistinguishable alternatives. When n = n' and m = m! there is 
only one alternative for photon absorption by the particles. In the case p = q, 
at difference with the rest of values p ^ q, there is not two-particle interference. 

We obtain the same conclusion, in a little bit more intuitive way, in the 
position representation. For instance, in the particular case I\l = Kr = AT, the 
joint detection probability becomes 

P{x , y) = P P ro(x , y) + \y K {x)\ 2 \v K {y)\ 2 cos ({q - p)(x - y)) (9) 

When the two momenta are equal the dependence on x~y 1 characteristic of the 
spatial two-particle interference, vanishes. 

This behavior can be understood in terms of the form of the initial state 
<po, which transforms for p = q into y/2\p >l \p >r that is a product state 
instead of an entangled one. Moreover, the normalization is incorrect: the 
normalized (in the single-mode sense) two-particle state is \p >l \p >r- We 
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have a discontinuity in the set of entangled states. For any p ^ q you have the 
entangled state <j> o. However, for p = q there is not an entangled state. 

The structure of the space of entangled states (with equal coefficients 1 /\/2) 
is not !ft 2 but !ft 2 — {p = q}- It is constituted by two disconnected parts separated 
by the line p = q. The same conclusion holds for the general states a\p >l |q >r 
+(3\q >l | P >r , where the mathematical space is S 2 (g) {5ft 2 — {p = g}} instead 
of S' 2 (g> 5ft 2 , with S 2 the Bloch sphere. 

This mathematical discontinuity translates into a physical discontinuity. If 
we evaluate the joint probability as the limit q —> p within the space of entangled 
states we have lim q ^, p P(x,y) = 2 P pro . This is the joint detection pattern for 
entangled momenta with q ss p (although p ^ q). It doubles at any point 
the pattern P pro of a product state. Consequently we have a discontinuity for 
equal values of the momenta. In its vicinity we have joint patterns with sharply 
different values. The transition between entangled and product states is not 
continuous. 


3 Multi-mode states 


We consider in this section the multi-mode case. In order to deal with analytical 
expressions we only consider Gaussian states, which have been extensively used 
in the literature (see [12] for non-Gaussian situations). The unnormalized initial 
state is 

$ 0 (p,g) = e - p2/Q2 e- q2/Q2 *e- p(l/p2 (10) 

The coefficients Q, Q* and P denote the spread of each exponential. 

When P~ 2 ^ 0 the above state is entangled. A good measure of the entan¬ 
glement degree for continuous variable systems is the Schmidt number [10, 11]. 
For our problem it can be easily evaluated analytically when 4P 4 > Q 2 Q 2 , 
giving S = (1 — Q 2 Q 2 /4P 4 ) -1 / 2 . The system is entangled when S > 1. 

The unnornralized wave function after the interaction, < F*(p, q), is evaluated 
using two times the Fourier transform (we omit some constant coefficients) 


®*(p,q) 



dye- ipx / h e- iq y/ R Y,Kbm 

n,m 


e i(p 0 +2nhK L )x/h e i(q 0 +2mhK R )y/h _ 


'y ' b n bm 


n,m 


dp 0 

dp 0 


dq 0 $ 0 (Po,qo) x 
dqo$o(po,qo) x 


S(p — po — 2 nhI<L)S(q - qo - 2ijiTiKr) = ^ b n b m $o(p ~ ZuUKl, q - ZriTiKr) (11) 

n,m 

The normalized wave function is $ = with N the normalization factor. It 

can be explicitly evaluated when 4P 4 > Q 2 Q 2 and is given by 

N~ 2 = J2 blbl l b r b s TTP 2 e-^ n2+r2)H2K ^ /Q2 e- i{m2+s2)h2K ^ /Q2 x 

n,m,r,s 
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e -4(mn+rs)h 2 K L K R /P 2 e a 2 Q 2 /& ex p 


/ 

V 


P 4 (/3 — fpr) 2 

2 ( 3 2 ~~ *) 


( 12 ) 


with a = ATiK L (n + r)Q 2 + 27 \K R (m + s)P 2 and /3 = AhK R (m. + s)Q* 2 + 
2TiK R {n + r)P~ 2 . 



Figure 2: P(o,q) in arbitrary units versus q in units Q = 1. The black, red and 
blue curves correspond respectively to P = oo, 1.1 and 0.75. 

We give next a graphical presentation of the detection patterns associated 
with these equations. We take as momentum unity Q = 1. In addition, we use 
the values Q* = 0.9, TiK l = 0.2 and hK r = 0.3. In order to see the variation 
of the patterns with the entanglement degree we consider three different values 
of P: oo,l.l and 0.75 that correspond to a product and two entangled initial 
states. We represent in Fig. 2 the two-particle detection pattern for p = 0, 
P(0, q) = ^(O, q) | 2 . In the evaluation we only consider the terms n = 0, ±1, ±2. 

We have that in the three cases the sum of multiple Gaussian terms (the 
terms in Eq.(ll)) reduces in a very good approximation to a single effective 
Gaussian distribution: P e ff{ 0, q) = a e ff exp(— q 2 /Q 2 e ff), with the values <r e // = 
0.7,0.65,0.42 for P = oo, 1.1, 0.75, and Q 2 e ff = 4000 for all the P's. Note that 
these single distributions are not normalized (the normalization is only for the 
sum over p of all the P(p, q)). 

Q e f f is a measure of the width of the distribution. More physically it quan¬ 
tifies to how many (p = 0,g)-modes effectively affect the diffraction. In all the 
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cases, for product states and for different degrees of entanglement, it is equal. 
On the other hand, cr e // represents the intensity of the diffraction process for 
the affected modes. In contrast with Qeff, it depends on the entanglement de¬ 
gree. In our example it decreases for increasing Schmidt’s numbers. For other 
values of p that behavior will be different (the total probability for all the p's is 
the same in the three cases). 

We have also studied the form of the two-particle diffraction patterns in the 
position representation. For the sake of shortness we do not present the complete 
analysis here. We only signal that, as in the momentum representation, a single 
Gaussian distribution fits well the multiple Gaussian terms sum. 


4 Identical particles 

In this section we study the multi-mode case when the two particles are identical. 
It is interesting because in addition to the interference and entanglement effects 
we must consider the exchange ones. The last effects are only relevant when 
the overlapping between the particles is large. Then the two particles must be 
diffracted by the same light grating (see Fig. 3). 



Figure 3: The same as Fig. 1 but with only one light grating. The mother 
particle decays into two daughter ones that interact with the grating. 

4.1 Entanglement and overlapping 

First of all we must analyze the relation between entanglement degree and over¬ 
lapping, the two measures that determine the intensity of the entanglement and 
exchange effects. The state of two identical particles (described by the entangled 
state <I>o previous to the consideration of the identity conditions) is obtained via 
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the standard procedure of (anti)synrmetrization: 

®ide(p,q) = {2±29)~ 1/2 {4> 0 { Pl q) ±4> 0 {q,p)) (13) 

with 9 the overlapping of the two particles 

f f ( 4P 4 - 0 2 0 2 \ 1//2 

d = J dP j dqn(p,q)Mq,p) = QQ* ( P 4 ( Q2 + Q 2) 2 Q , Qt ) ( 14 ) 

The last expression is valid for P 4 > Q 4 Qi/{Q 2 + Q 2 ) 2 (for instance, for 
the values used in the previous section). When this condition does not hold 
the overlapping cannot be evaluated analytically. The signs + and — in ± refer 
respectively to bosons and fermions. 

Through this paper we have not considered the spin (or electronic) part of 
the quantum state because it is not relevant in the diffraction dynamics. When 
the spin variables are taken into account the (anti)symmetrization refers to the 
complete state and, for example, the spatial wave function of fermions can be 
symmetrized. By assuming that the two particles are in a symmetrized spin 
state, |s >i |s >2 or (|s >i | — s > 2 +| — s >i |s > 2 )/-\/2, we do not need to 
consider these cases here. With this choice the spatial part of the state must 
be symmetrized for bosons and antisymmetrized for fermions. The extension to 
antisymmetrized spin states is immediate. 

Although there is some controversy on the characterization of entanglement 
in systems of identical particles, following the criterion in [13] it is clear that $ide 
is entangled. For instance, for fermions, the state is non-entangled if and only if 
4>id e is obtained by antisymmetrizing a factorized state [13]. Thus, for fermions 
(13) represents an entangled state. Moreover, in this approach one clearly sees 
that entanglement is associated with <f>o, the state that is (anti)symmetrized. 
Then it is natural to assume that the amount of entanglement in is the 
same of $o- As signaled in the previous section we use the Schmidt number to 
quantify the entanglement contained in <f>o- 

Two interesting conclusions can be easily obtained: 

(i) We fix Q and Q * and vary P over all the values showing entanglement, 

{QQ*/ 2) 1 / 2 < P < oo. When P —> (<2Q*/2) 1//2 we have 9 —> 0 and S —> oo. The 
overlapping is very small, whereas the entanglement tends to very high values. 
On the other hand, for P —>• oo these limits are 9 2QQ*/{Q 2 + Ql) and 

S —> 1. It is simple to see that this value of the overlapping is maximum. Now 
the system tends to a factorized one whereas the overlapping is finite, reaching 
the maximum value compatible with Q and Q*. We observe a complementary 
behavior. The entanglement is maximum when the overlapping is null and vice 
versa. 

(ii) Symmetric entangled states represent a special case in the relation be¬ 
tween entanglement and overlapping. For symmetric entanglement, Q = Q*, 
the overlapping reaches its maximum value 9 = 1 independently of P. This is 
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so for any type of symmetric entangled state, & sy m{p, q) = ^sym(qiP), because 
f dp J dq\Q sy m{p, q)\ 2 = 1 is due to the normalization. For symmetric entangled 
states we have maximum overlapping regardless of the particular entanglement 
degree of the state. 

In addition, symmetric states have another relevant property. When we take 
<E>o symmetric is undefined (it has the form 0/0) for fermions. No pair of 
fermions can be prepared in a symmetric entangled state (with our choice of 
the spin part), a property that can be seen as a natural extension of Pauli’s 
exclusion principle from product to non-factorizable states. In effect, when the 
state is a product one, <f>(p)<f>*(q), the symmetry condition implies <j> = (/*, a 
relation that is forbidden by Pauli’s exclusion principle. 

On the other hand, for bosons, we have = < f>g ym (p, q) that corresponds 
to a state without symmetrization (in the sense of symmetrization for identical 
particles). The exchange effects vanish: in the presence of symmetric entangle¬ 
ment there are not exchange effects in two-boson systems. This result confirms a 
previous analysis in [14], where a similar behavior in the position representation 
was found for pairs of non-entangled identical bosons in a two-slit arrange¬ 
ment: in the cases where the overlapping between the two bosons is large (it 
is increased by the diffraction process to values close to unity) the two-boson 
diffraction pattern is indistinguishable from that of a product state and does 
not show the typical characteristics associated with exchange effects. 

4.2 Diffraction patterns 

After analyzing the relation between S and 9 in the initial state we evaluate the 
diffraction patterns. By similitude with Eq. (11) the final state is 

$ ide(p,q ) = bnbmi&oip — 2nhK , q — 2mhK) ± <j>o(q — 2nhK,p — 2m%K)) 

n,m 

(15) 

The calculation of the normalization factor M is a little bit more involved than 
in the previous section. In 1| 2 we have four terms, two direct and two crossed 
ones. The integration of each one of them over the full momentum space gives a 
contribution to the squared normalization factor. Moreover, we must take into 
account that the crossed terms must be added to the direct ones in the case of 
bosons, but subtracted in that of fermions. For instance, the contribution of 
the first term is 

•AC 2 = E KKJbrb.i rC 
n,m,r,s 

e -4(m 2 + s 2 )K 2 K 2 /Ql e -4,(mn+rs)h 2 K 2 /P 2 e y. 2 e/4 exp 



C _ j e ~4(n 2 +r 2 )h 2 K 2 /Q 2 
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where £ -2 = Q~ 2 + Q~ 2 , n = 2TiK(2mQ~ 2 + 2rQ~ 2 + (s + n)P~ 2 ) and JI = 
2hK(2sQ~ 2 + 2nQ~ 2 + (m + r)P~ 2 ). This expression is valid when P 4 > £ 4 . 
For other values of the parameters an analytical expression cannot be obtained. 
It is simple to see that the two direct (D) terms are equal, and the crossed (C) 
ones also. Finally we have M 2 = 2Mp ± 2A Iq. Md is given by N .* and Me by 
a similar expression with obvious modifications in the order of the coefficients 
n, m, r, s. 



Figure 4: The same as Fig. 2 with identical particles. The black, red and blue 
lines correspond to the cases P = 200,1.1 and 0.75, whereas continuous and 
dashed ones represent respectively bosons and fermions. 

Next, we represent these patterns. With the same values for Q and Q* of Fig. 
2, K = 0.2 and for the parameter P the values 200 (negligible entanglement), 
0.75 (small overlapping) and 1.1 (intermediate case) (see Fig. 4). We observe a 
clear difference between bosons and fermions. In the first case the distribution 
is peaked around the point (0, 0), whereas in the second this point is forbidden 
by Pauli’s exclusion principle. Instead, we see a valley around this point with 
two peaks around. In addition to this general trend we see that the visibility 
(difference between maximum and minimum values) is strongly dependent on 
the entanglement degree of each curve. 
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5 Conclusions 


We have shown that entangled states are diffracted by a light grating in a differ¬ 
ent way than product ones. This is the essence of multi-particle interference. In 
an alternative view, our results reflect the dependence of the light-matter inter¬ 
action on the separability of the state. In our proposal one can easily visualize 
how the joint detection patterns change with the entanglement degree. 

The discontinuity found in the single-mode approximation provides, in prin¬ 
ciple, a method to experimentally study the disconnected character of the space 
of entangled states. It would be also interesting to test the cancellation of ex¬ 
change effects by symmetric entanglement described above. This result, which 
corroborates previous findings for non-entangled states [14], suggests the ex¬ 
istence of a rich and unexplored relation between entanglement and particle 
identity. For instance, the complementarity between overlapping and Schmidt’s 
number reflects a type of competition between both effects. 

In this paper we have restricted our considerations to thin light gratings. 
For thick ones we move to the Bragg regime of the Kapitza-Dirac effect. One 
expects to find similar effects to those described here, but this must be verified 
by an explicit calculation. 

We have not discussed the experimental realization of the proposal. As the 
diffraction regime of the Kapitza-Dirac effect has already been tested [5, 6], we 
must concentrate on the sources of entangled particles for the arrangement. For 
two gratings, a decaying particle initially at rest, seems to be a good candidate. 
The single-mode case can be approximated by restricting the two generated 
beams to some particular directions using collimators. For only one grating, the 
photodissociation of molecules traveling towards the standing light wave, looks 
to be a more adequate process than the decay of unstable particles. 

In the case of electrons a nanotip source can be interesting. The antisym- 
metrization properties of such a source have been studied in [15]. They lead to 
a dip in the joint detection probability. Although the state of electrons appears 
to be entangled in momentum (Eqs. (15) and (16) in [15]), it is not true entan¬ 
glement: the state is obtained by (anti)symmetrization of two-particle product 
states (Eqs. (12) and (13) in [15]) and then, according to the criterion in [13], 
the entanglement is associated with the (anti)symmetrization procedure. How¬ 
ever, due to the Coulomb interaction between the electrons, there can be true 
entanglement in the system. An extended analysis including that interaction 
would be necessary to evaluate the entanglement and to see the viability of an 
experimental test. 

Acknowledgments I acknowledge the support from Spanish Ministerio de 
Ciencia e Innovation through the research project FIS2009-09522. 
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